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Abstract 

The scaling limit as T ^ of the antiferromagnetic three-state Potts model on the 
square lattice is described by the sine-Gordon quantum field theory at a specific 
value of the coupling. We show that the correspondence follows unambigously from 
an analysis of the sine-Gordon operator space based on locality, and that the scalar 
operators carrying solitonic charge play an essential role in the description of the 
lattice model. We then evaluate the correlation functions within the form factor ap- 
proach and give a number of universal predictions that can be checked in numerical 
simulations. 
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The study of antiferromagnetic models is a notoriously difficult problem of statistical 
mechanics. The dependence on the lattice structure produces a variety of behaviours 
much richer than in the ferromagnetic case and forces a case-by-case analysis. Many 
antiferromagnetic models possess a critical point (often at zero temperature), so that 
their continuum limit can also be investigated through field theoretical methods. Among 
them, the square lattice three-state Potts model has been the object of both numerical 
and theoretical studies. While this model has been known for longtime to be critical 
at zero temperature, the issue of the approach to criticality presents several subtleties, 
including the identification of the correct scaling variable. Very recently, the authors of 
reference exploited a mapping onto a height model to identify the excitations on the 
lattice at nonzero temperature and explain the anomalous corrections to scaling previously 
observed in Monte Carlo studies 0. 

It is the purpose of this note to point out that the scaling limit of the square lattice 
three-state Potts antiferromagnet is exactly solvable due to its equivalence to a specific 
point of the sine-Gordon model (an integrable quantum field theory), and to derive from 
this continuum approach some universal predictions that can be checked through simula- 
tions on the lattice. 

As it often happens for two-dimensional systems exhibiting a gaussian critical point, 
a relation of the scaling limit with the sine-Gordon model is expected. The actual task is 
that of understanding what the symmetries and the operators of the lattice model become 
in the field theoretic language. Here we will show that these identifications follow in a 
quite natural way from the analysis of the sine-Gordon operator space based on locality, 
and that the topologically charged scalar operators play a major role in the continuum 
description of the lattice model. Matrix elements for these operators were ffist computed 
in 0, where they arise in the description of the scaling Ashkin- Teller model. We will 
then illustrate how the correspondence at the operatorial level translates in the language 
of particle excitations, in order to exploit the exact S'-matrix solution of the sine-Gordon 
model and extract from that the correlation functions for the operators of physical interest. 
In particular, we will give a number of specific and universal predictions suitable for a 
numerical verification on the lattice. 

The three-state Potts model is defined by the lattice Hamiltonian 

H = -J^5,^^s,, (1) 

where the spin variable Sj on the j-th lattice site takes the values (colours) 0,1,2 and 



1 



the sum is over nearest neighbours. The Hamihonian is invariant under permutations of 
the colours. Here we are interested in the antiferromagnetic case (J < 0) on the square 
lattice. 

At zero temperature, the system is in one of the ground states in which nearest neigh- 
bours spins have different values. This problem is equivalent to the three colouring prob- 
lem of the square lattice which admits an exact mapping onto a specific point of the 
six- vertex model 0, ^. From the exact solution of the latter we know that our zero- 
temperature system is critical and that its long-distance behaviour is described by a free 
massless boson. Relevant operators that can be identified through the zero-temperature 
analysis on the lattice ^, ^, ^ are the staggered magnetisation = (— 1)-J'i+j''2 e2i7rsj/3 
with scaling dimension 1/6, the uniform magnetisation aj = e^*'^'^^/^ with scaling di- 
mension 2/3, and the staggered polarisation^ = (— 1)-'^^"'^ Z]i(25si,s^ — 1) with scaling 
dimension 3/2. We identify the j-th site of the square lattice through a pair of integers 
(ji;j2), and call even (odd) sublattice the collection of the sites with ji + j2 even (odd). 

At non-zero temperature, the model develops a finite correlation length and its scaling 
limit is described by the perturbation of the gaussian fixed point through the continuum 
version of the thermal operator Sj = J2i^si,Sj- This perturbed theory is the sine-Gordon 
model defined by the euclidean action^ 

= J <9aV99"v9 — /iCOS/3(y9^ , (2) 

for some value of the coupling (3 to be determined. The theory (H) describes the scaling 
limit of several lattice models whose critical point corresponds to a conformal field theory 
with central charge C = 1 (see e.g. Refs. 0, |10| for other examples discussed in the 



framework of this note). In order to proceed with the description of the antiferromagnetic 
Potts model we need to recall few points about the operator content of the sine-Gordon 
model. 

At criticality (yU = 0), the boson field can be decomposed into its holomorphic and 
antiholomorphic parts as (p{x) = (f){z)+(f){z), where we introduced the complex coordinates 
z = Xi + ix2 and z = Xi—ix2- The scaling operators of the theory are the vertex operators 

Vp^p{x) = e*W(-)+P'^(-")l , (3) 



^The primed sum indicates summation over the next nearest neighbours of j. 

^Throughout this note we use the standard notation f3 for the sine-Gordon dimcnsionless couphng. 
No confusion with the inverse temperature of the lattice model should be made. 
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with conformal dimensions (A, A) = {p'^ /87T,p'^ /8tt) and spin s = A — A. They satisfy 
the gaussian operator product expansion 

We see from this relation that taking around Vp^^p^^O) by sending z —>■ ze^'^'^ and 

z —>■ ze~'^^'^ produces a phase factor e^*'^'''^'^, where 

71,2 = ^ {PiP2 - P1P2) (5) 

is called index of mutual locality. If 71,2 is an integer the correlator {Vp^^p^{x)Vp^^p^(fi)) is 
single valued and the two operators are said to be mutually local. Since 71,1 = 2s, the 
operators which are local with respect to themselves (the only ones we are interested in 
here) must have integer or half integer spin. 

In the off-critical theory (Q), the operators wich are local with respect to the perturbing 
operator cos [3ip ~ Vp^p + V^p^^p form a 'local sector' into which all the operators of interest 
for the description of the lattice model are expected to fall. This locality requirement 
selects the operators Vp^p with p — p = Anm/ [3, m integer, namely 

Vp{x) = Vp,p{x) = e*^"^(^) , (6) 

and 

Un,m{^) = V.,,..^.p 2.^(x) = e4^'^(")+T-'^(")] , n = 2s = 0,±l,..., m = ±l,... 

2m /3 2m /3 

(7) 

Here we introduced the 'dual' boson field which is (p{z) — (f){z) at criticality and satisfy 
the relation 

dxa "'^dXjS ^ ^ 

The operators Vp and f/o,m are scalars (s = 0) and have scaling dimensions Xp = p^ /An 
and Xo^m = vrm^//5^, respectively (X = A + A). The action (0) describes a relevant 
perturbation of the gaussian fixed point for < Svr; in this range the only operators f/o,m 
which are relevant (Xo,m < 2) are those with |m| < 3. 

The lowest operators with \s\ = 1/2, i.e. \E' = f/±i_i and \E'* = U±i^-i, are complex 
conjugated two-component spinors with conformal dimensions A„ m given by 



8 1^47r 
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A±i,i = A±i,_i = - ( ^ ± 2 + — , (9) 



^±i.i = ^±W=o (l^T2 + ^ . (10) 
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It is known since Coleman [0] that the sine-Gordon model is equivalent to the theory 
of a Dirac fermion with four fermion interaction, the Thirring model. The expression for 
concides with the bosonisation formula for the Thirring fermion originally derived by 
Mandelstam |12|. The value (3"^ = An for which the dimensions (P), (0) take the free 
fermionic values corresponds to the free point of the Thirring model. 

It follows from these considerations that the integer m in (0) is a fermionic charge. 
Since the Thirring fermions correspond to the solitons interpolating between adjacent 
vacua of the periodic bosonic potential, we will also call m topologic charge. Hence, we 
can make a distinction between neutral scalar operators Vp{x) with a non-zero vacuum 
expectation value, and charged scalar operators Uo^m{x). 

Going back to the Potts model it is not difficult to identify the continuum limit of 
the lattice operators among the sine-Gordon operators. Both the staggered and uniform 
magnetisation are charged with respect to colour permutations and must be found among 
the f/o,m- Comparison with the known scaling dimensions gives ~ ^^0,1(2;) and 

(r{x) ~ f/o -2(3^)5 and selects 

/3 = , (11) 

as the value for which the sine-Gordon model describes the scaling limit of the three- 
state Potts antiferromagnet on the square lattice. This immediately fixes the scaling 
dimension of the thermal operator S{x) ~ cos (3 if (x) to be 3/2, in agreement with the 
result obtained in [jT| by studing the vortex excitations on the lattice. 

On the basis of these identifications we can see how the symmetries of the lattice 
model translate in the sine-Gordon language. The group 6*3 of colour permutations can 
be decomposed into the Z3 transformations associated to cyclic permutations plus the 
complex conjugation of e^*'^^j/^. The latter operation simply corresponds to the complex 
conjugation of the sine-Gordon exponentials, while the elementary Z3 transformation 
maps into the shift (p ^ + 2ti / (3. The charge coincides with the topologic charge m 
(mod 3). 

The lattice operators are also characterised by their parity under the transformation 
which exchanges the even and odd sublattices: £j and aj are even, while Ej and Vj are 
odd. It appears that in the continuum limit this parity property corresponds to (— 1)™. 
The staggered polarisation Vj is invariant under colour permutations and must correspond 
to the operator ?7o,3 + f/o,-3! which has indeed the expected scaling dimension 3/2. We 
summarise the situation in Table 1. 
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Table 1 



The equivalence with a particular case of the sine-Gordon model ensures that the 
scaling limit of the lattice model is an integrable quantum field theory whose associated 
scattering theory is known exactly [|T3|. The elementary excitations are a pair of con- 



jugated particles and A_ (the sine-Gordon soliton and antisoliton) . The fact that 
= Gtt falls in the repulsive sine-Gordon regime > Air ensures that no other particles 
are present in the spectrum. It can be interesting to mention that the Za-preserving fu- 
sion A^A+ A_ is forbidden here because it violates the sublattice parity (—1)™; it is 
instead characteristic of the scattering theory of the ferromagnetic case [|1^, in which 
the lattice plays no role. 

Due to the factorisation of multiparticle scattering in the integrable quantum field 
theories, the scattering theory is completely determined by the two-particle S'-matrix 
defined by the relation^ 

Aai9i)A,{e2) = E &1 - 92)A,{92)A,{9i) , a,b = ±. (12) 

c,d=± 

The non-zero scattering amplitudes are given byQ |13[ 

SlX{9) = SZZ{9) = Soi9) , (13) 

sinh^ 

SUB) = Szm = ^-^^^^SM, (14) 

sinh ^ 

= = (15) 



The rapidity 9 parameterises the on-sheh energy and momentum of a particle as {p^,p^) = 
(m cosh 0, m sinh 9) . 

^For reasons that will become clear later, it is useful to give the results referring to generic values of 
P in the sine-Gordon model, being understood that the scaling limit we are dealing with corresponds to 
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with 



/•oo sinh f ( 1 — ^ ) Qj- 

So{e) = -exp{-t — ^ ^ sin — ^ , (16) 

Jo X sinh|^cosh| vr 



From the S'-matrix one can compute the form factor^ 



f = ^- (17) 



„.(»!,. ..,e„) = (0|t(0)|A„,(ft),. ..,/!„.(»„)), a, = ± (18) 

which in turn determine the spectral decomposition of the correlation functions. The 
form factors satisfy the equations ^7\, |18[ 

^*,...,a.,ai+i,...,a„(^l) • • • , 6**, ^'i+l, . . . , 9^) = 

F*_,J^i + 2m,^2,...,en) =e^^'^^*-F*_^^^^^ , (20) 

where 7$,a in the last equation is the index of mutual locality between the scalar operator 
$(x) and the soliton (for a = +) or the antisoliton (for a = —). Since these particles are 
created by the operators Uo^±i{x), equation (^ gives for the scalar operators 

P 

7yp,± = lVp;Uo,±i = ±^ , (21) 
7c/o.™;± = 1Uo,m;Uo,±i = • (22) 

Equations (0) and (pO]) can be used to determine the 'lowest' form factors for the 



operator $, i.e. the non-vanishing matrix elements ([18|) with the smallest n (n > 0) 
which fix the initial conditions of the form factor bootstrap^. The lowest form factors for 
the operators of interest in this note are 

pnsh ^ 
smh ^(6^12 - «7r) 



^-^2n),...,-.,M(^i,---,^H) =CHn^o(%) , (24) 



-sg{m),...,-sg{m)\ 

i<j 



^|0) is the vacuum state. 

^The higher form factors contain additional pairs A^{6)A_{d') in the asymptotic state and are related 
to the lowest ones by a recursive equation associated to particle- antiparticle annihilation [|l7t. 
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where 9ij = 6i — 6j, sg{m) denotes the sign of m, and the Cm are normahsation constants. 
The function 



e \ r^dx sinh [| (l - J) 

— t Sllili '^^^■^ 

satisfies the equations 



• 2 {iiT—d)x 

Fo{9) = -I sinh - exp M — ^ ^"^ ) (25) 

2 ^\Jo X sinh g cosh f sinhx ^ ^ ' 



Fo{9) = Soi9)Fo{-e) , (26) 
Foie + 2t'K) = Fo{-e) . (27) 

The lowest form factors determine the first term in the large distance expansion of the 
correlation functions. Using the subscript c to denote the connected correlators we have 

{V,{x)VM)c = /^^|F]:^_(0,,0,)|2e-^N(-^h^^+™^t^^^) 



27r 271 

+ 0(e-^^'l"l) , (28) 

ml J 271 ZTT 

+ 0(e-(™+')^''l"l) , (29) 

where M is the mass of the soliton and we took m > 0. 

There are a number of simple and distinguished universal predictions of the sine- 
Gordon description that can be tested in numerical simulations of the square lattice 
three-state Potts antiferromagnet. The 'exponential' and 'second moment' correlation 
lengths C,^ and ^l"'^ associated to an operator $(a;) are defined as 

(<l>(x)$*(0))e ~ exp(-|a;|/e$) , |a:| ^ oo , (30) 

,2nd ^ A Id'x\x\^^x)^*m^ 

^* 1^4 jd^xmx)^*io)), ) ■ ^ ^ 

The leading behaviour of these quantities nearby the critical point is 

(32) 

(33) 

where = 1/(2 — X^) = 2, and t measures the deviation from criticality. It appears from 
numerical simulations [§] and has been discussed in that for this zero-temperature 
critical point the scaling variable is t = e"^/'^^. Our operator identifications and equations 
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, (ID then give ^^ = l/M and 

U/h = I, (34) 

fs/h = I, (35) 

fv/h = I, (36) 
/|"7/e ^ 1, (37) 
/fV/, ^ 0.439. (38) 

The last two ratios are evaluated in the lowest form factor approximation which is known 
to give estremely accurate results for integrated correlators^. We estimate that our error 
on these two quantities does not exceed 1%. 



The universal ratios (|3^)-(p8D have not yet been measured in simulations. Lattice 
estimates of these quantities would provide the direct confirmation that the universality 
class of this antiferromagnet is described by the sine-Gordon field theory with /? = a/Gtt 
and the operator identifications of Table 1. In this respect, it is important to mention 
that the study of the scaling limit of the model on the lattice is complicated in practice 
by strong corrections to scaling ||l|, Cardy et al. argued in [Q] that such corrections 
are originated by a strictly marginal operator which couples to the temperature giving an 
effective thermal dependence to the stiffness K of the bosonic field, which is related to 
the sine-Gordon coupling as (3 = Q^/K. Their estimate of this (nonuniversal) dependence 
for Kef fit) corresponds to I3lff{t) = 67r - 20.9(4)t - 37.8(2)^^ + . . .. Hence, according 
to this analysis, the leading corrections to scaling can be fitted by letting the critical 
exponents and the critical amplitudes vary with j3eff- Concerning the lattice verification 
of our predictions, the ratios (|51D-(P7D are not affected by this effect since their values in 
the sine-Gordon model are determined by symmetry properties of the operators which do 
not depend on the coupling^. The ratio (/^"°'//(T)e// is instead an increasing function of 
(3effi SO that the asymptotic value (|38D should be approached from below as t — > 0. For 
example, we find 0.38 when computing this ratio at Plffit = 0.1) ~ Gtt — 2.47. 

Acknowledgments: I thank John Cardy for helpful discussions and remarks. 
Note added: When this work had been completed reference appeared which is 



''See e.g. |l^ for similar computations in other sine-Gordon related statistical models. 
*At least in the range 0^ > 4tt in which there are no bound states. 
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devoted to the form factors of the sine-Gordon topologically charged operators and whose 
main purpose is the determination of their "conformal" normahsation. 
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